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ABSTRACT 

We study the effect of introducing a weak antiferromagnetic interplanar ex- 
change coupling in the two dimensional frustrated Heisenberg model. We show 
that a ferromagnetic(FM) ordering of chirality - i.e., same chirality on adjacent 
planes - is energetically favoured, thus leading to bulk violation of the discrete 
symmetries parity(P) and time reversal(T). 
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Ever since the discovery of high T c superconductors [1] and the observation of 
their layered nature, it has been conjectured [2 — 4] that these materials are de- 
scribed by a ground state that explicitly violates the discrete symmetries parity(P) 
and time-reversal(T) macroscopically. Theoretical interest in these ideas began 
with the work of Kalmeyer and Laughlin [4], who (approximately) mapped the 
Heisenberg model on a triangular lattice, to a bosonic FQHE problem at filling 
fraction v = 1/2, with semionic excitations. An apparently very different line 
of investigation was initiated by Affleck-Marston [5] and Kotliar [6], who intro- 
duced the notion of 'flux phases'. Working with a frustrated Heisenberg model 
on a square lattice, Wen, Wilczek and Zee [7] generalised the 'half-flux' phase of 
Affleck-Marston to a 'quarter-flux' phase which they called the chiral spin liquid 
(CSL). They found that this state explicitly violated P and T macroscopically and 
in the low energy long wavelength limit, its effective action led to semionic statis- 
tics, thus corroborating the Kalmeyer-Laughlin picture. More recently, Laughlin 
and Zou [8] have shown that the Gutzwiller projected CSL state is identical to 
the Kalmeyer-Laughlin state, paving the way to a three dimensional generalisation 
of the physics underlying the FQHE. The concept of flux phases has also been 
extended to the doped situation and generalised flux phases have been shown to 
be plausible ground states of the doped t-J model [9]. On the experimental front, 
many novel experiments were both suggested [10] and performed [11] to look for P 
and T violation, which appeared to be a robust prediction of all anyonic theories. 
But the experimental situation remains confused in the face of conflicting evidence. 

Much of the earlier theoretical work was confined to studies of single planes. 
But lately, there have been several attempts to extend flux phase ideas to the 
fully three dimensional situation [12]. However, it is also of both theoretical and 
experimental importance to incorporate weak three dimensionality -i.e. , to study 
the effect of weak interlayer couplings -in planar phenomena. We focus on this 
particular aspect in this letter. We study the effect of a weak antiferromagnetic 
interlayer spin-spin coupling (well motivated by neutron scattering studies [13]) on 
two dimensional CSL ground states. By perturbatively computing the correction 
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to the ground state energy, we show that a FM ordering of chirality on adjacent 
planes is preferred. Our work is close in spirit, but somewhat complementary to 
the work of Rojo and Canright [14], who studied the ordering of anyons on adjacent 
planes when a static scalar potential is introduced between them. But our work 
studies the ordering of the ground states of a microscopic model -the frustrated 
Heisenberg model- whereas the starting point of their work involves a gas of anyons. 
The two calculations, therefore, cannot be directly compared, since it is not yet 
possible to explicitly derive a gas of anyons from the frustrated Heisenberg model 
or any other microscopic model. 

Let us consider 2p planes, each of which has a spin S = 1/2 sitting on the 
sites of a square lattice, with Heisenberg antiferromagnetic interactions between 
all nearest neighbours within each plane. The ground state of this model is well 
known to be Neel ordered. However, one of the important effects of doping this 
model with mobile holes is to induce frustrating interactions [15], which cause an 
instability towards generalised flux phases. Qualitative features of these phases are 
captured by the CSL states of the frustrated Heisenberg (J-J') model given by 

2p 2p 

Ho = JJ2 E S?-S? + J']T £ S?-S? (1) 

a=l <i,j>£n.n a=l <i,j>£n.n.n 

which is also an interesting model in its own right. We shall use this model as our 
starting point. In Eq.(l), i is the two dimensional site index common to all the 
planes and the index a identifies each plane. As argued in Ref. [7], the CSL state, 
characterised by the order parameter < Sj • Sj x >, where i,j and k are the 
vertices of an elementary triangle, is a local minimum of this model for sufficiently 
large J'. In fact for slightly modified Hamiltonians, it is a plausible ground state. 
It is this ground state which has anyonic excitations and motivates the study of a 
gas of anyons which forms the basis of theories of anyon superconductivity [16]. In 
this letter, we shall focus our attention on the mean field description of this CSL 
state. 
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We use a fermionic description for the spins given by 

Sr = £4k/^ and Sf a = J2di* a pd^ (2) 

a,j3 a,/3 

where we have distinguished alternate planes by the nomenclature of the fermions 
as c and d planes. Also, for clarity and brevity of notation, we shall henceforth drop 
the index a and the summation over a which now goes from 1 to p, since odd and 
even planes have been distinguished. In terms of the fermions, the Hamiltonian Hq 
can be rewritten, after a Fierz transformation followed by a Hubbard-Stratanovich 
transformation as 

H 0= \Xij c ia c 3<* + xfjAjja + h.C.] + ^[ a 0i(4a c ^ ~ 1) + 4(4a^a ~ X )] 

{i,j},a i,a 

+ 1 v y c ' d K c4 + — v Y c,d V' d 

<i,j>£n.n <i,j>£n.n.n 

(3) 

where Xif are the Hubbard-Stratanovich fields. The notation {i,j} in the first 
summation stands for summation over both nearest and next nearest neighbours. 

Following WWZ [7], we introduce the mean field ansatz for the chiral spin 
liquid state for both the c and d planes. For the nearest neighbour links, 

<X-;L>=^ /4 and (X&j) = ge-**'* (4) 

where i here, is a site on the odd sublattice. However, for the diagonal links, the 
WWZ ansatz admits a two-fold degeneracy corresponding to the two possible chi- 
ralities - i.e., the flux through each elementary plaquette, which is now a triangle, 
could either be positive or negative. Since we wish to study the ordering of chi- 
ralities on different planes, we allow for independent chiralities on the two planes. 
Thus, the n.n.n links are described by 

(Xi,i-x+y) = 5if c ' d and (Xi,i-x-y) = -kf A (5) 

where Si = +(— )1 for % belonging to the even (odd) sublattice. f c = f d implies 
that the chiralities on the adjacent c and d planes are the same (FM ordering) and 
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f c = —f d implies an AFM ordering. In the absence of any interplanar coupling, the 
two possibilities obviously remain degenerate. Finally, for the Lagrange multiplier 
fields we have 

(<4> = (4> = o, (6) 

so that the fermions are no longer subject to the 'no double occupancy' constraint 
at each site. The constraint is now enforced only on the average. 

Notice that the mean field ansatz for the CSL state divides the square lat- 
tice on each plane into two sublattices. Thus, we may take the spatial Fourier 
transformations separately for the odd and even sublattices, with respect to a 2-d 
wave vector k which now runs over the reduced Brillouin zone (RBZ). Hence, the 
momentum space mean field Hamiltonian is given by 

H MF = E ^ih^L + <lhUL] (7) 

k€RBZ,a 

where tp^a = ( C ka> C D and ^ka = (^ka'^LJ' and e stand for odd and even 
respectively and 

k - U 

2f c ' d [cos(k x + k y ) — cos(k x — k y )] 2gr[ — i cos(h x ) + cos{k y )} 

2g[i cos(k x ) + cos(k y )] —2f c,d [cos(k x + k y ) — cos(k x — k y )} 

(8 

This Hamiltonian can be diagonalised by a unitary transformation yielding 

tdiae = ( ,,,/) 



E k 



with E^ d = (|Ak| 2 + {^]f) 2 ) 1 ^ 2 in terms of the transformed variables (7k a ?7ka) 
and (r]^ a , rj£ ) for the valence band(V) and conduction band(C) fermions in the c 
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and d planes respectively. The ground state has the valence band completely filled 
in both the planes and its energy, in terms of the mean field variables, is given by 

£» MF = § E E 9 2 + }E E z 2 - 2 E ( 10 > 

ieodd j(i)£n.n ieodd j(i)£n.n.n k€RBZ 

In the absence of any interplanar coupling, the FM and AFM orderings of 
chirality remain degenerate. To lift the degeneracy, we introduce a weak interlayer 
Heisenberg antiferromagnetic coupling given by 

H int = J"J2^^t (11) 

i 

where Sf and Sf refer to the spins on the c and d planes respectively. Such an 
interaction is particularly appropriate for the copper oxide systems and leads to 3-d 
Neel ordering in the undoped insulating phase. J" has been estimated from neutron 
scattering experiments [13], to be about five orders of magnitude less than the in- 
plane coupling J. We treat Hi nt as a static perturbing potential between the two 
species of fermions on adjacent planes. This is accomplished by taking momentum 
space Fourier transformations with respect to a 2-d wave vector. Thus, despite the 
extension of the problem into the third dimension, inter-layer particle transfers are 
avoided and the essential layered nature of the original problem is retained. The 
relative weakness of J" with respect to J justifies this approach. 

In the fermionic representation, 

J" t + 

Hint = — C ia C iP d lp d ia ( i2 ) 
i,a,/3 

which, when Fourier transformed with respect to 2-d wave vectors, becomes 

Hi n t = — ^ [ C k+qa C k/3 rf J-q/^k'a + C k+qa C k/3 rf Xq/3 rf k'a] • ( 13 ) 
k,k',q€i?BZ 

Notice that a change in momentum in the c plane is compensated by an opposite 
change in momentum in the d plane. We now evaluate the total ground state 
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energy, treating Hi n t as a perturbation, for the two cases of FM and AFM orderings 
of chirality. The unperturbed ground state is given by 



(Ground state) = ]J^\0) ® ( 14 ) 

and the unperturbed ground state energy is given in Eq.(lO) . The FM ground 
state has f c = f = f d , whereas the AFM ground state has f c = f = —f d , so 

that the Boguliobov transformation coefficients and hence the definition of the 

V c 

transformed fermions r] k ' a differ in the two cases. It is now straightforward to 
rewrite Hint in terms of the transformed fermions and compute Epm — Eafm- 

At first order, we find that 

k 

since e^/E^ is odd under reflection about the k y -&xis and the summation over k 
includes both positive and negative k x . This result is easily understood, since at 
first order, the only term in Hint that contributes involves no momentum transfer q 
hence, the two planes are essentially independent and the degeneracy between FM 
and AFM orderings of chirality is not lifted. At second order too the degeneracy 
is not lifted. We find that 

p (2) _ p (2) /j^_\2 ST \ , Ck'-q e k ' w ek+q 6 k ^ 

^FM AFM ^1AT> L-i TP TP ^ TP TP ^ TP TP > ^ TP TP >' 

Zl\ ^ i^k + -frk' + -^k+q + -&k'-q -&k'-q -&k' -^k+q ^k 

(16) 

Making the changes k — > k — q, q — > — q, and k' — > — k' succesively in the dummy 
variables, and using e_k = e k , A_ k = A k and E^ = E-^, we find that 

E FM ~ E AFM = ~( E FM ~ E AFm) = °- ( 17 ) 

However, the third order contribution does lift the degeneracy and is given by 

E FM ~ E AFM = E A + E B + E C (18) 
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where 



J" 3 ^ , ek+q + ek+q+q' ~ &\l ~ ekCk+qek+q+q' 

ek' + ek'+q' - e^-q - ek'ek'+qek'-q 

Ek + -Ek+q+q' + Ek'+q' + Ek'-q 
_ e k'-q + ek'-q-q' ~ ~ ek'ek'-q-q'Ck'-q j 
4(£'k + -Ek+q+q' + -Ek' + Ek'-q-q') 



£ J" 3 1 

B ~ 2N 3 ^ / (E k + E k+q + E k ' + Ek'-q) (E k + E k+q+ q' + Ek'+q' + Ek'-q) 

[ek+qek'(^k 5 k+q+q'^k'+q'^k'-q + h.C.) 

+ ekek'-q(5k+q^k+q+q'5k'5k'+q' + h.C.) 

+ ek+q+q'ek'+q'(5k^k+q^k'^k'-q + h.C.)], 



(20) 



and 



Ec = Il y i 

l, v i + ^ k +q + + E k '- q )(E k + E k+q+q ' + E k ' + Ek'-q-q') 

[ek+qek'-q(5k^k+q+q' 5 k'5k'-q-q' + 
+ ekek'(^k+q^k+q+q'^k'-q^k'-q-q' + ^- C 
+ ek+q+q'ek'-q-q'(^k+q^k' 5 k'-q + h.C.)]. 

(21) 

(Here, ek = ek/Ek and 5^ = Ak/Ek for any momentum k.) The k- summations in 
Eqs.(19), (20) and (21) were performed numerically using a Monte Carlo routine, 
for different values of J'/ J, with the corresponding values of / / J and g/J being 
obtained by minimising Eq in Eq.(10) with respect to / and g. Our numerical 
results are tabulated below. 
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Table I 



V 1 T 
J / J 


tl J 


„ 1 T 

91 J 


1 77>(3) ip(3) \ 1 j 
\ FM nj AFM)/ J 


0.50 


0.015 


0.23 


- 1.7 xlQ- 4 N(J"/J) 3 


0.55 


0.024 


0.24 


- 2.9 xlQ- 4 N(J"/jf 


0.60 


0.035 


0.24 


- 4.6 xlO- 4 N(J"/J) 3 


0.65 


0.046 


0.23 


- 6.8 xlO- 4 N(J"/J) 3 


0.70 


0.060 


0.23 


- 8.2 xlQ- 4 N(J"/J) 3 



Thus, for any value of J'/J for which the CSL state is a local minimum, and 
for J" > 0, (which is the case for copper oxides), the FM ordering of chirality is en- 
ergetically favoured. Notice that the energy difference is an extensive quantity and 
scales linearly with N. In fact, using typical values for La2CuO^ ( J = 1200°^ 
and J" = 0.03°K), and assuming iV ~ 10 15 , Eafm - E FM ranges between 3.2°fT 
and 15.4°K for J' j J between 0.5 and 0.7. Thus, despite the weakness of the inter- 
layer coupling, its potency is effectively increased by its extensivity. Hence, at low 
enough temperatures, the weak Heisenberg antiferromagnetic interlayer coupling 
has sufficient strength to tilt the scales in favour of a FM ordering of chirality. 

In our calculation, we have completely ignored gauge field fluctuations -i.e., the 
phase fluctuations of the order parameter Xij an d the fluctuations of the Lagrange 
multiplier field aoi- These fluctuations could lead to a substantial contribution 
to the ground state energy. However, they cannot lift the degeneracy between 
the FM and AFM orderings of chirality, since they only act within each plane. 
Thus, as long as these fluctuations do not destabilise the mean field ground state, 
Efm — Eafm and consequently, the ordering of chirality is determined only by 
the interplanar coupling. 
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We have worked within the framework of the J-J' model, which is a limiting 
case of the t-t'-J-J' model. However, we expect the qualitative aspects of our re- 
sult - i.e., the tendency towards FM ordering of chirality - to be valid even for the 
generalised flux phases of the doped t-J model, at least for low doping. Notice that 
our result suggests that bulk P and T violation is an inescapable consequence of 
CSL ground states of models that are relevant to high T c superconductors. More- 
over, despite the controversy regarding the observation of local P and T violation, 
bulk P and T violation has certainly been ruled out in the cuprate compounds [11]. 
Hence, our calculation disfavours models with CSL ground states as candidates to 
describe the doped Mott insulating phases of the copper oxides. 
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